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In quenched QCD, where the dynamic creation of quark-antiquarlc pairs out of the vacuum is 
neglected, a confined baryonic system composed of three static quarks exhibits stringlike behavior 
at large interquark separation, with the formation of flux tubes characterized by the geometry of 
the so-called Y ansatz. We study the fluctuations of the junction of the three flux tubes, assuming 
the dynamics to be governed by an effective bosonic string model. We show that the asymptotic 
behavior of the effective width of the junction grows logarithmically with the distance between 
the sources, with the coefficient depending on the number of joining strings, on the dimension of 
spacetime and on the string tension. 



PACS numbers: 12.38.Aw, ll.25.-w, 21.30.Fe 

I. INTRODUCTION 

One of the crucial features characterizing quantum 
chromodynamics (QCD) is confinement: the fundamen- 
tal constituents of strongly interacting matter (quarks 
and gluons) are not observed as asymptotic states, and 
the low-energy hadronic spectrum consists of colorless 
states only. 

Because of its nonperturbative nature, a formal proof 
of confinement from first principles has so far been elu- 
sive, and even the degrees of freedom responsible for this 
phenomenon are subject to debate. However, some low- 
energy properties of hadrons and QCD forces — see, e.g. 
Ref. [1] for a review — can be accurately modeled in 
terms of an effective string picture [2-6] , which describes 
the infrared properties of hadrons in terms of a fluctu- 
ating, thin (almost unidimensional) flux tube joining the 
color sources. At sufficiently large interquark separations 
the lowest-energy excitations of the confined system are 
associated with collective degrees of freedom correspond- 
ing to transverse stringlike vibrations of the flux tube, 
whereas the excitation spectrum of the gauge degrees of 
freedom inside the tube is much higher lying. 

While it is unproven that low-energy aspects of con- 
fining quantum field theories can indeed be explained in 
terms of an effective string theory with universal features, 
strong evidence in favor of this conjecture is provided by 
recent comparisons between results from lattice simula- 
tions and bosonic string predictions (surveyed below). 
The observation of a linearly rising potential between 
color sources separated at distances R does not automat- 
ically imply the existence of an effective string descrip- 
tion, however, the universality of the subleading 1/i?- 
coefficient found in several gauge models [6-18] provides 
such a nontrivial test. At present not many quantitative 
predictions that exceed the classical limit are available 
to compare lattice data to. In this article we increase 
the number of such nontrivial predictions by calculating 
the width of "baryonic" flux-tube junctions for general 
geometries. 

Note that while we label our configurations as bary- 



onic, we expect our predictions only to apply to pure 
Yang-Mills theories with static external charges. The 
bosonic string model is unlikely to be a good approxima- 
tion to the baryons of real QCD with sea quarks, which 
are likely to decay into baryon-meson pairs, before the 
string limit of large distances can be reached. 

For the simplest physical scenario, i.e. a ("mesonlike") 
pair of static, infinitely heavy, confined color sources, 
the effective model has been developed since the early 
1980s [2-5] . Later on this description was reformulated in 
terms of an expansion about the long-string vacuum [19]. 
More recent theoretical developments include Refs. [20- 
22], and are reviewed in Ref. [23]. The large-distance 
string behavior has been observed in numerical lattice 
simulations of the torelon spectrum (corresponding to 
closed strings) and static potentials (corresponding to 
mesonic open strings) of SU(3) lattice gauge theory [6-9] 
as well as of various other gauge models [10-18]. 

The effective string picture also predicts the width of 
the flux tube to grow logarithmically as a function of the 
interquark distance as well as the coefficient of the log- 
arithm [4]; this was addressed and confirmed in various 
numerical lattice studies [15, 24-30]. 

Lattice simulations of the baryonic setup [1, 31-34] in- 
dicate that the flux-tube profile interpolates between the 
so-called A geometry at short distances (where the effec- 
tive one-gluon exchange dominates), and the Y ansatz 
for separations between the sources of the order of or 
larger than approximately 0.8 fm. This Y ansatz that is 
relevant in the infrared region is characterized by a junc- 
tion where the flux tubes meet (see Fig. 1). The corre- 
sponding leading order string corrections to the baryonic 
potential have been worked out in Ref. [35]. 

Building upon the procedure used in this reference, in 
this article we study the width of the junction, assuming 
that this is generated by string fluctuations of an effective 
string theory with the lowest dimensional term given by 
the Nambu-Goto (NG) action [36, 37], with string tension 
a. We consider the leading order nontrivial behavior in 
the limit of large separations between the static color 
sources. 



FIG. 1: A snapshot of the fluctuating flux tubes (ideahzed as strings) of three static color sources Q, joined at a common 
junction; ip denotes the position of this junction, relative to its classical location minimizing the total string length. The 
components of parallel and normal to the plane containing the static color sources are also displayed. 



II. CALCULATION SETUP 



In formulas: 



The calculation is performed in £)-dimcnsional Eu- 
clidean spacetime with D — \ spatial dimensions of infi- 
nite extent and a periodic time coordinate t £ [0, T). We 
consider the general case of 2 < n < Z? static "quarks" 
spanning a (n — l)-dimensional hyperplane. We assume 
the action to be minimal when the n strings meet in a 
common junction. For n = 3 this is indeed the case, 
unless one of the angles of the triangle defined by the 
three sources exceeds the critical value 2tt/3. In this lat- 
ter case the junction will be fixed to the position of the 
corresponding source and the system can be decomposed 
into two mcsonic strings. Note that also for many n > 3 
geometries the classical configuration will be character- 
ized by different geometries, with two or more distinct 
junctions (see also Ref. [38]). In principle our calculation 
can be extended to these cases. 

During their time evolution, the strings span n dif- 
ferent world sheets (see Fig. 2); each of these blades is 
bounded by the (straight) worldline of a static quark on 
one side, and a generic worldline spanned by the fluctu- 
ating junction on the other side. 

Classically, the ground state fulfills the constraint of 
the minimal area of the string world sheets. Therefore 
the position of the junction is determined by the require- 
ment of minimal total string length. Furthermore, the 
balance of tensions means that a X]a=i = which 
implies equal angles between the strings at any time, 
denote unit vectors along the direction between the junc- 
tion and the quark a. Assuming that the string dynamics 
is described by the NG action S'ng means that the quan- 
tum weight of a generic configuration of a string world 
sheet is proportional to the exponential of its total area. 



5, 



NG 



(2.1) 



where the position on the surface of the string world 
sheet is parametrized by the coordinates C,i,i € {1,2}. 
Although this bosonic string model is nonrenormalizable 
(because it is nonpolynomial) and anomalous (except in 
26 spacetime dimensions) [5, 39, 40], it can be consid- 
ered a legitimate starting point in the construction of an 
effective theory, and it can be shown to agree with the 
Polchinski-Strominger effective model [19] up to and in- 
cluding the next-to- leading order [21]. Let X'^((^), with 
fi € {1, . . . , D}, denote a map from the world sheet to 
the spacetime, embedding the world sheet; the induced 
metric of Eq. (2.1) is given by 
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(2.2) 



To proceed with the quantum treatment, it is natural 
to fix the reparametrization and Weyl invariance to the 
"physical" gauge, allowing us to describe the transverse 
displacements (for our purposes, the Weyl anomaly can 
be neglected, because it vanishes in the limit of large 
distances). This means that only transverse fiuctua- 
tions ^a(t, s) of the string world sheets a E {l,...,n} 
around the classical configuration are considered as phys- 
ical. The time t and parameter s label the position on 
string world sheet (blade) a. In particular, s denotes the 
spatial distance from the worldline of the quark a, i.e. the 
classical position of the junction is given by s = L^. The 
junction worldline ip{t) fiuctuates within the hyperplane 
spanned by the quarks (changing the minimal area of the 
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FIG. 2: Left-hand side: World sheets spanned by the fluctuating strings during their time evolution. Right-hand side: Surface 
of one of the string world sheets. 



blades) as well as in the D ~n remaining transverse spa- 
tial directions. From continuity, we obtain the boundary 
conditions for the transverse fluctuations ^a{tys)' 

^a{t,La + ea-V>{t)) = V±a{t), (2.3) 

where, (p±a = '^ — s-ai^a'V)- The transverse fluctuations 
^a{t, s) vanish at the location of the quarks (s = 0), and 
are periodic in the time t, with period T. 

For technical reasons we assume the junction itself to 
have a finite mass m. This results in a static energy and 
in a kinetic term. The parameter value m should not af- 
fect the large-distance results, La(T 3> m, that we present 
below, and indeed it cancels in the calculation. Expand- 
ing the NG action around the equilibrium configuration 
yields 

(2.4) 

where again (^i , C2 are world sheet parameters and 

5|| = (T ^ (^LaT + jdtea- v(t)^ = ctLyT . (2.5) 

Ly = J2a ^bovc dcuotcs the total string length. (Note 
that Ba = 0.) In the T — > cxd limit the string thick- 
ness can be extracted from the partition function, which 
to leading nontrivial order is given by 

Z = e-(-^-+")^ J exp (-f /^i f[ Za{^) , 

(2.6) 



where Za{}p) denotes the partition function for the fluc- 
tuations of a given blade that is bounded by the junction 
worldline tp(t): 

Z,(v) = y 2)^, exp (^-| 1 \dia\^ . (2.7) 

The string partition functions Zaij^) arc Gaussian fimc- 
tional integrals and can be calculated as follows: 

Zj<p) =e-f/l^€™-lVet(-Arjr(^-2'/\ (2.8) 

where ^min.a is the minimal-area solution for given ^'(t). 
Ar„ denotes the Laplacian acting on the domain (blade) 
To- $min,a(ij'S) is harmonic and satisfies the boundary 
conditions Eq. (2.3) [35]. Below we will evaluate this ex- 
pression to the leading order in terms of the fluctuations 

In contrast to the mcsonic setup, the world sheets F^ 
are in general no rectangles. However, the determinant 
in Eq. (2.8) can still be calculated by decomposing the 
boundary tp(<) of F^ into a sum over Fourier modes and 
conformally mapping the resulting domains to rectangles, 
as shown in Appendix A. Carrying out this mapping and 
taking the limit T — > cxd, Jahn and de Forcrand [35] de- 
rived the subleading term of the [n = 3) baryonic poten- 
tial Vqq,: 
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a 



1/^ 



. a<b 



(2.9) 



r 



We confirm this result. For tiie equilateral case L = Li 
L2 = Eq. (2.9) simplifies to: 



V, 



qqq,A 



(L) = 3aL - 



D-3tt 



such that the approximation coth(C) 1 holds. Neglect- 
ing subleading terms in L, the second part of integration 
gives 
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y + 0(L-2). (2.10) 



dw 



1 



(J ji^ mup' + naw coth(i/;L) 



— InL , 



III. STRING THICKNESS AT THE JUNCTION Therefore, to leading order, the result reads 



The bosonic string model yields a prediction for the 
thickness of the fluctuating strings. The width of the 
junction itself can be calculated by taking the expecta- 
tion value: 



(3.1) 



The action S is defined in Eq. (2.4) and can also be 
read off from the partition function Eq. (2.6). We split 
the string width into contributions (y^^), perpendicular 
to the hypcrplane spanned by the n quarks, and ((^"^) 
within the hypcrplane of the quarks: 



(3.2) 



Note that ip^ lives within a (D— n)-dimensional subspace 
and yj" fluctuates in the remaining n — 1 spatial direc- 
tions. This differs from the definitions with respect to a 
given blade a of Sec. II above, (p±a [{D — 2)-dimcnsional 
fluctuations] and ip^^a (one-dimensional fluctuations). 

In the T — > cx) limit, the perpendicular contribution 
reads 



1 r°° 

71" Jo 



1 



Ha COth(wLa) 



(3.3) 

We present more details of the calculation in Appendix B. 
For n strings of identical length L = Li = L2 = ■ ■ ■ = Ln 
we obtain 



1 f° 
(D-n)- / 

Jo 



dw 



mw"^ -t- naw coth(wL) 



(3.4) 



We split this integral at w = C/L. where C is an arbi- 
trary L- independent constant. The first part of integra- 
tion is subleading in L. We can choose C large enough. 
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D-n 1 L 

In — 

n TTcr Lq 



(3.5) 



(3.6) 



where we have absorbed an arbitrary constant into Lg; 
the width of the junction, orthogonal to the plane 
spanned by the quarks, grows logarithmically with the 
distance. Lo will depend on Z?, n, and the microscopic 
(ultraviolet) details of the gauge model. 

We can perform a consistency check by comparing the 
above result to the mesonic case (n = 2). The logarith- 
mic behavior of the width of a flux tube connecting two 
quarks in the string picture was predicted many years 
ago by Liischer, Miinstcr and Weisz [4]. The result they 
found for Z? = 4 is 
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7ri\/2 A 



(3.7) 



where denotes the string tension and A represents a 
cutoff scale. The effective string width was studied again 
by Caselle et al. [26] and by Gliozzi [41], calculating the 
deviation of the transverse coordinates of the string from 
the respective Green function. The result they obtained 
for the mean squared width Wq in D = 3, determined at 
the symmetry point of the string world sheet is 



27rcr Rc 



(3.8) 



where a denotes the string tension, R the interquark dis- 
tance and Rc is an ultraviolet scale. 

We divide the string connecting quark and antiquark 
into two parts of equal length (up to small longitudinal 
fluctuations), connected in the middle by a junction. We 
can then apply Eq. (3.6) for £» = 3 and £» = 4. The 
above predictions indeed coincide with our results where 
we identify a ^ hP, L^2L' ^ 2R and Lq = 2A = 2i?c- 

Now let us turn to the string width within the plane 
of the quarks. This is only well defined for n > 3 and the 
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most interesting case is the n = 3 baryon. In contrast to 
the perpendicular width (cp-^^), our calculation applies 
to n = 3 only since we assume the sources to lie in a two- 
dimensional plane. In this case the n > 3 minimal string 
configuration will usually contain more than one junc- 
tion, unless junctions are fixed at the positions of quarks 
and do not fluctuate. It turns out that our calculation of 
((pll^) cannot easily be generalized to higher dimensional 
planes, i.e. to ri > 3. In the baryonic equilateral case we 
obtain [Eq. (B23)]: 



\^qqq.A/ SanJo rh / L + {w - a / L^) coth{w) ' 

(3.9) 

where m = 2m/{3a) and a = {D - 2)/(127r(T). The 
result for nonequilateral configurations can be obtained 
from Eq. (B19) below. 

We split the integral in analogy to the above discussion 
of the perpendicular fiuctuations. However, one finds 
that, like in the calculation of the baryonic potential, a 
pole emerges. If one views the NG action as the first 
term within an effective string theory then this pole has 
to be canceled by counterterms arising from the inclusion 
of higher dimensional operators. In this sense it should 
not affect the leading order result. Assuming this, the 
parallel contribution to the width of the junction turns 
out to be: 

where again subleading contributions are suppressed and 
Lc is an undetermined constant. Note that the coefficient 
above is by a factor 4/(Z? — 3) larger than the one in 
front of the logarithm within the expression for (tp"*"^) of 
Eq. (3.6). One factor 2/{D — 3) corresponds to the ratio 
of independently fiuctuating parallel over perpendicular 
components while another factor of 2 is expected from the 
stronger restoring force for perpendicular displacements, 
relative to parallel ones. 

IV. CONCLUSIONS 



with a Steiner junction can be obtained from Eq. (3.3). 
The width within the plane of the sources also grows log- 
arithmically as a function of the separation and we have 
calculated this for n = 3. The result for the equilateral 
case is displayed in Eq. (3.10) while the general result 
can be calculated from Eq. (B19). We also confirm the 
result of Ref. [35] for the baryonic potential, Eqs. (2.9) - 
(2.10). 

The masonic flux-tube width has already been inves- 
tigated in lattice simulations of different gauge theo- 
ries [15, 24-30]. While most of these studies confirm 
the logarithmic broadening of the mesonic string, with- 
out much lattice spacing dependence, it should be noted 
that Ref. [25] found such broadening only at fixed lattice 
spacings but the string width to actually shrink, possi- 
bly to zero, if the continuum limit was taken. This is 
also incompatible with the result of Ref. [42] of a vanish- 
ing overlap between a thin string state and the ground 
state wave function. Further lattice studies are required 
to resolve this controversy. 

The question if and at what distances our string pre- 
dictions become valid can be addressed by lattice simu- 
lations of baryonic configurations in SU(3) gauge theory 
at large L in D = 3 and D = A spacetime dimensions. 
While this is numerically quite challenging, at least the 
simplified case oi D ~ 3 Z3 gauge theory can be mapped 
to a two-dimensional Potts model, allowing for precise 
numerical simulations [34, 43]. These show consistency 
with the potential of Eq. (2.9). 



Acknowledgments 

We thank Oliver Jahn for providing us with the de- 
tailed notes of his calculation of the baryonic potential. 
We also thank O. Jahn and Philippe de Forcrand for en- 
lightening discussions. M.P. gratefully acknowledges fi- 
nancial support from the Alexander von Humboldt Foun- 
dation and from INFN. The University of Regensburg 
hosts the Collaborative Research Center SFB/TR 55 
"Hadronenphysik mit Gitter-QCD" . 



In this article, we studied the width of the junction of 
flux tubes in baryonlike systems composed of infinitely 
heavy, static color sources (quarks) at large distances L 
from the junction. Assuming the low-energy aspects to 
be governed by the dynamics of the bosonic Nambu-Goto 
string model, we have shown that the width of the junc- 
tion grows logarithmically with the distance between the 
quarks. In particular the quadratic width orthogonal to 
the {n — l)-dimensional plane spanned by n equidistant 
quarks [the baryons of SU(n) gauge theories] in D > n 
spacetime dimensions reads [Eq. (3.6)]: 

{V^') = In^. 4.1 

This also applies to (and generalizes) the mesonic case 
(n = 2). The corresponding result for general geometries 



APPENDIX A: CONFORMAL MAPPING OF A 
WORLD SHEET TO A RECTANGLE 

Here we provide the ingredients for the calculation of 
the fiuctuations at the junction. We follow Ref. [35], con- 
formally mapping the blade (see Fig. 2) to a rectangle. 
The minimal-area solution for a fixed position of the junc- 
tion, ^min,ait, s), is harmonic and satisfies the boundary 
conditions Eq. (2.3): 

A|,„i„.a = , ^min,a{t, La + Ba ■ (fit)) = ¥'_La(0 ■ 

The determinant in Eq. (2.8) is computed with Dirichlet 
boundary conditions on the domain = {(t, s)|0 < s < 
La ~t~ Ba ■ ip(t)}. In terms of the Fourier components 
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of (p{t), |,nin,a IS given by 



1 



sinh(ws) 



VT ^ smh(wLa) 

where w = 2'Kn/T. The integral in Eq. (2.8), which 
represents the change in the minimal area due to the 
transverse fluctuations tp^a, can now be calculated: 

.2. V- di 

mm, a ^smin.a j_i / t \\ 1 2 

'^^A^~Q^, '^F~ " 2^WC0t\\{wLa)\ip^^A.a\ 



(A3) 



The determinant in Eq. (2.8) is obtained by mapping the 
domain Fa conformally to a rectangle L'^ x T. Note that 
the conformal map fa{z) ^ z + J^w'^waS^'^ has to be 
complex differentiablc. Its coefficients c^a ai'c fixed by 
the constraints: 

faiK+lt) = La + Ba ■ ifit) + it + 0{ip^) . (A4) 

One easily sees that L'^ ~ La + -^^a ■ Vo- To leading 
order in the conformal map is then given by 

.fa{z) - - + ^ E -^TTT-,^'"' + ^(^') • (A5) 

Vr ^„ smh(u;La) 

This conformal mapping changes the Laplacian by a 
scalar factor: 



At 



e^^-WAi^xT, Pa{z) = -lln\d,fa\^ . (A6) 



The variation of the determinant of the Laplacian with 
respect to a holomorphic mapping of T onto some other 
region F via the function f{z) can be calculated by means 
of the Alvarez-Polyakov formula (see e.g. Ref. [2]): 



In 



det(-Ar) 1 



det(-Ap) 127rjgr 
+ ^ J^dhdAn\dJ\^d, ln|9,/|2. (A7) 

Here z(t) is an arbitrary parametrization of dV and z' = 
dz/dr. In our case the first integral above vanishes and 
thus, from the conformal map Eq. (A5), we obtain to 
leading order 

/ d^zdjn\dja\^d,ln\dja\^ 

= E'«'|ea-<^^|'cotli(wLa) + 0(^'), (A8) 

w 

where we used the fact that the Fourier coefficients sat- 
isfy (fi-n, = if'^. In Ref. [2] the rectangle with peri- 
odic boundary conditions in time is further conformally 
mapped onto a circle resulting in: 



(A9) 



where 7y(r) denotes the Dedckind 77 function. Collecting 
the above results, we obtain for the determinant of the 
Laplacian with respect to the blade a: 

iT 



(AlO) 



det(-ArJ = r;^( — 



X exp 



coth{wLa)\ea ■ V^pj 



APPENDIX B: CALCULATION OF THE WIDTH 
OF THE JUNCTION {(p^) 

In this appendix we provide more steps for the calcu- 
lation of Eqs. (3.3) and (3.9). The calculation is carried 
out for n quarks located in a plane. This configuration 
with only one junction might not be stable for more than 
three quarks [38]. However, the result for the orthogonal 
contribution can easily be generalized to configurations 
of n quarks, distributed in a (n — l)-dimensional hyper- 
plane. 

The thickness of the string at the junction can be cal- 
culated taking the expectation value of [g^c Eq. (2.4)]: 



To do this, we have to consider integrals 
Dip exp 



(Bl) 



a— 1 \ i 



D 



247r 



coth(wLQ)|ea • ifv 



(B2) 



We can replace the integral in the first term above by a 
sum over Fourier components: 



dt\(p\ 



(B3) 



We denote the plane that is spanned by the spatial unit 
vectors of the n strings = {ea,xi Sa,y, 0, . . .) as the x-y 
plane. These n unit vectors obey the relation = 0. 

The x-y components of ip (or any other vectors) carry the 
superscript "||". We obtain 

^coth{wLa)\ea ■<pi\'^ 

a 

= \vi,xf^el^xCoth{wLa) + |(/7i^j^|2Eea,yCoth(u;La) 

a a 

+2 (Re(v^ll , JRe((^lL,,) + Im(^![,JIm(^ll ,,y)) 

X E ^".^^".y coth(wLa) . (B4) 
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We define , Ay and as 



Thus we obtain 



, a {D- 2)w^ 

A^ = [ -w+ — 

2 247r 

A„ = — H 

,'(T (L'-2)ii;3 



XI ^1^^ COth(ti>La) 

CI 



(B5) 



^ea.xea,yCOth{wLa) 



Note that A^+Ay as well as the combination A^Ay — 
are invariant under rotations within the x-y plane. Since 
the angles between the n strings are equal, we can 
parametrize the unit vectors by 

Ba = (cos(27ra/n), sin(27ra/n), 0, . . .) . (B6) 



Ax Ay A^^ — 



(7 D-2W3y . / (a-b)\ 
+ - — — > sinM27r-!^ ^ 



X coth(z«La) coth(wL6) . 



(B7) 



In the case of ri — 3, the geometrical coefficients are 
ctab ~ 3/4. This results in the simplification: A^Ay — 



3 / a 



^^247!"" ) J2a<b'^0th{wLa)c0th(wLb). 



Al 



We can split the integral over in Eq. (B2) using 
Iv'u'.J-aP = IVii'P ^ \Vw ■ Gap into parts that are parallel 
and perpendicular to the plane of the quarks: 



Dip exp 

Tlip^ exp 
X / D<^" exp 



n — ~ 2 ^ — A Q 

+ coth(wLa)|ea • ip^ 



It! \ a / 



+Im(¥'ll,i.)Im(V5i,y)) ^Ro 



Xcoth(i«La) + \^UAx + I^IL + 2 (Re(^ll ,jRe(^ll 

a J 

(B8) 

I 



Here, are the 13 — 3 components of tp^ that are per- where 
pendicular to the plane spanned by the quarks. We ab- 
breviate the first functional integral above as Ii and the „ f 
second as l2- The solutions of these Gaussian integrals ^J- = / Tup-^ip-^"^ exp < — - > | mw"" + aw 
read J I 



in- 



D-3 



h 



h = 



n 



+ aw coth{wLa) 



n 



u>>0 

where we defined 



^1^2 - 4AL 



, (B9) 
(BIO) 



In = /Dv?IIv'II exp 



mw + (TU) 



X coth( WLa 

a 

(B13) 

COth('inLa^ 



Xcoth(u;La) , 



Al — Cyj — 2A^ , 

A2 = Cm — 2y4.y . 

We are interested in the expectation value 



(Bll) 



(B12) 



+2(Re{ipi^jReiipij+lmi^ijImiiplj)A^^] I .(B14) 



Let us recall that cp-^ is {D — 3)-dimcnsional. Performing 
the first integral yields 



\'Uj'>0 / tu>0 
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so that 



h T 

^ ■u)>0 



(D-3) 



With w = 27rn/T, in the Umit of large T we obtain 
TT Jq mw"^ + (Jw coth(wLa) 



(B16) 
1 

(B17) 



is calculated analogously. We perform the integral /|| 
[Eq. (B14)]. In terms of and A2 defined in Eq. (Ell), 
one is left with 



^ n 



E ■ (BI8) 



\u)>0 / t(;>0 



The parallel contribution to the width of the junction Therefore, combining /j| with Eq. (BIO) yields 
I 



h 2 



A1+A2 4 1 



- — / dw- 



1 



77iit; + (1 — auP')Ci 



^^AiA2-AAl^ 3 an Jo w m'^w^ + 2mw{l - aw^)Ci - Aaw^ {Cif + + aw^)^C2 ' 

M 



where 



2m _D'~2 
3 cr ' 127r(7 ' 



(B20) 



(B21) 



Ci = i ^coth(wLa) , 

a 

C2 = ^^aabCoth{wLa)coth{wLb) . (B22) 



For an equilateral baryon we have X)a<6 '^^b = | ^-nd 
L = Li — L2 — L3. Thus, for this special case, we 
obtain 



dw 



Scttt io w2m/L+(u;-w3a/L2)coth(w) ■ 

(B23) 



a<b 
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